Robust half-metallicity and metamagnetism in Fe z Coi_ x S2 
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The Fe^Coi-z^ system is predicted, on the basis of density functional calculations, to be a 
half metal for a large range of concentrations. Unlike most known half metals, the half metallicity 
in this system should be very stable with respect to crystallographic disorder and other types of 
defects. The endmember of the series, C0S2, is not a half metal, but exhibits interesting and 
unusual magnetic properties which can, however, be reasonably well understood within the density 
functional theory, particularly with the help of the extended Stoner model. Calculations suggest 
strong electron-phonon and electron-magnon coupling in the system, and probably a bad metal 
behavior at high temperatures. 
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God made the integers, all else is the work of man 

(L. Kronecker) . 

Half metals (HM), materials that are metals in one spin 
channel and insulators in the other, are attracting sub- 
stantial interest recently, mostly because of potential ap- 
plication in spintronics devices, e.g., spin valves jL|. Un- 
fortunately, although a few dozen various materials have 
been predicted to be HM on the basis of band structure 
calculations, there are hardly any that have been convinc- 
ingly confirmed to be such by an experiment. A possible 
exception is Cr02 where a spin polarization of up to 90% 
was measured by the Andreev reflection technique 0] (in 
some other materials there is indirect evidence such as 
integer magnetic moment or optical spectra consistent 
with half-metallic bands). The usual explanation of such 
a discrepancy between the theory and the experiment is 
that half-metallicity in these materials is very sensitive 
to crystallographic disorder and stoichiometry (see, e.g., 
Ref. f||). Indeed, no materials have been predicted to be 
HM in a wide range of concentrations of the constituents, 
and insensitive to disorder. 

In this Letter I point out one such material. Namely, 
I show that the pyrite alloys Fei-^Coa^ are HM for the 
most of the concentration range (0.1< x < 0.9), and not 
sensitive to ordering of Fe and Co in the metal sublat- 
tice. On the other hand, at x > 1 |Q, calculations predict 
magnetic collapse under pressure, and a metamagnetic 
behavior just before the collapse. Both half metallicity 
and metamagnetism can be explained by competition be- 
tween the kinetic (band) energy, and the Stoner (Hund) 
interaction. 

Experimentally, the system of pyrite solid solutions, 
(Fe,Co,Ni,Cu)(S,Se)2, is amazingly rich. FeS2 is a non- 
magnetic semiconductor, in agreement with the band 
structure calculations HH]. With as little as 0.1% Co, 
by some data 0, it becomes a ferromagnetic metal, 
and remains ferromagnetic all the way through C0S2, 
and further on until approximately Coo.9Nio.1S2. It was 
noted that the magnetic moment, M, per Co atom in 
Fei_ 2 ,Co 2 ,S2 solid solutions stays close to 1 /is in an ex- 
tremely wide range from x « 0.1 — 0.2 to x = 0.9 — 0.95. 
To the best of my knowledge, no explanation of this fact 



has been suggested till now. At larger x, M decreases 
to w 0.85 [Ib @]- Starting from C0S2, one can also sub- 
stitute S with Se. The Curie temperature, Tq, rapidly 
decreases, and magnetism disappears at Se concentra- 
tion of 10-12%, while M decreases only slightly ||. At 
larger Se concentrations, the material shows metamag- 
netic behavior with no sizeable spontaneous magnetiza- 
tion, but with magnetization of 0.82—0.85 /ib appear- 
ing abruptly when applied magnetic field in exceeds ap- 
proximately (220a;5 e — 25) tesla. A very similar behavior 
was observed in pure C0S2 under pressure || , suggesting 
that the magnetic effect of Se is just the density of states 
(DOS) reduction. 

As shown below, all these effects are reproduced by 
the standard local spin density (LSD) calculations, and 
find their explanation within the extended Stoner model 
(ESM) ||. Further substitution of Co by Ni leads to a 
Mott-Hubbard transition into an antiferromagnetic state, 
which is, in contrast to the system considered here, 
poorly described by the LSD. Further doping by Cu leads 
to a superconductivity, presumably due to high-energy 
sulfur vibrons |^p| . 

To understand the behavior of the Fei_ x Co x S2 sys- 
tem, I performed several series of density functional LSD 
calculations: First, virtual crystal approximation (VCA) 
was used in conjunction with the Linear Muffin Tin Or- 
bital (LMTO) method jyj. Then, I did several calcu- 
lations using rhombohcdral supercells of 4 or 8 formula 
units. Finally, I checked the results against more accu- 
rate full-potential linear augmented plane wave calcula- 
tions for pure C0S2, FeCoaSg, and (Feo.25Coo. 75)82 
in VCA. The results were consistent with the LMTO cal- 
culations. All calculations were performed in the experi- 
mental crystal structure of C0S2. In reality, the S-S bond 
in FeS2 is 4% longer. The effect on the band structure is 
not negligible, particularly near the bottom of the con- 
ductivity band |l3|] . However, the difference is not impor- 
tant for the purpose of the current Letter, namely the half 
metallicity of Fei_ 2: Co a ;S2 alloys, and understanding the 
basic physics of its magnetic phase diagram. The result- 
ing magnetic moments are shown in Fig|l]. In good agree- 
ment with the experiment, the magnetic moment per Co 
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is exactly 1 hb for the Co concentrations 0.3< x < 0.9. 
The same holds for the rhombohedral supercell calcula- 
tions (Fig.|l]), demonstrating stability of the HM state 
with respect to crystallographic disorder. I performed 
calculations for ordered FeyCoSig, FeaCoSg, FeCoS4, and 
FeCosSs, and found that already FeyCoSig (x — 0.125) 
has magnetic moment of 1 /is per Co. Although the orig- 
inal paper Q implied that the total magnetic moment 
resides on Co, this is not true: for instance, in FeyCoSig 
less than 30% of the total magnetization (0.45 /j,b) re- 
sides on Co. The nearest neighbor Fe (6 per cell) carry 
w 0.15 fin each. 8% of the total moment resides on S, 
about the same relative amount as in C0S2 [Q. So, one 
cannot view the low-doping Fei_ a; Co a ;S2 alloys as mag- 
netic Co ions embedded in polarizable FeS2 background, 
as for instance Fe in Pd. 

The best way to understand the physics of this alloy is 
to start with FeS2- FeS2 is a nonmagnetic semiconductor 
with a gap between the tig and e g states (^|| . The reason 
for that is that sulfur forms S2 dimers with the ppa states 
split into a bonding and an antibonding state. The latter 
is slightly above the Fe e g states and thus empty. The 
other 5 S p states are below the Fe ti g states, hence the 
occupation of the Fe d bands is 6, just enough to fill the 
narrow t2 g band. Magnetizing F2S2 would require trans- 
fer of electrons from the t2 g into the e g band, at an energy 
cost of the band gap A per electron. The Stoner param- 
eter /, which characterizes the gain of Hund energy per 
one electron transferred from the spin-minority into the 
spin-majority band, appears to be smaller than A « 0.75 
eV (For the conductivity band in FeS2, due to hybridiza- 
tion with S, / is smaller that in pure Fe |l5| and is 
0.55 eV). However, if we populate the same band struc- 
ture with x <C 1 additional electrons per formula unit, 
we can either distribute them equally between the two 
spin subband, or place x' in the spin majority band and 
x" — x — x' in the spin majority band (the total magnetic 
moment in \xb is then M = x' — x"). In the latter case 
we gain a Hund energy of — JM 2 /4, but we lose kinetic 
(band structure) energy because some states have only 
single occupancy and thus more high-energy states have 
to be occupied. Using effective mass approximation for 
the conductivity band one can express this kinetic energy 
loss in terms of concentration x and magnetic moment 
per Co /3 = M/x as A{x/2f' 3 [(l + (3f/ 3 + (1 - /3) 5/3 ], 
where A — 3ti . (-Arr) 2 ^ m* is the effective mass, and 

10m* V 4-7T V ' 



V is the volume per formula unit. 

Minimizing the total energy with respect to f3 gives the 
equilibrium magnetic moment in ESM: 

13 = (5A/3/)[(l + P) 2 ' 3 - (1 - p) 2 / 3 ]/(4x)^ 3 

The solution is a universal function ^(Ix 1 ^ 3 /A). Note 
that (3(z) = for z < 10 • 2 1 / 3 /9, and P{z) = 1 for 
z > 5/3. Correspondingly, with doping the material re- 
mains paramagnetic until x reaches x\ — 2(10^4/9/) 3 , 



and then the magnetic moment per Co atom gradually 
grows until the concentration reaches the X2 = (5A/3I) 3 . 
At larger dopings the material remains half-metallic with 
/3 = 1. Eventually, DOS starts to deviate from the ef- 
fective mass model and the polarization /3 may become 
smaller than 1 again. One can estimate the critical 
concentrations X1.2 using the following values, extracted 
from the band structure calculations for FeS2: m* w 0.8, 
/ « 0.55 eV. This gives x\ ss 0.16 and X2 ~ 0.26, in qual- 
itative agreement with the experiment and LSDA calcu- 
lations (Fig.|l|). 
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FIG. 1. Experimental and calculated magnetic mo- 
ment per Co in Fei-^Co^^ alloys. VCA: virtual crystal 
approximation, s/cell: supercell calculations, experimen- 
tal data are from Ref. M. 



This explains the on the first glance unexpected re- 
sult: in contrast to most known half metals Fei_ a; Co 2 ;S2 
remains a HM for a large range of concentrations, is in- 
sensitive to crystallographic disorder, and probably not 
very sensitive to the state of the surface either: the be- 
havior qualitatively described by the universal function 
(3 above is determined by the competition of two large 
energies: The band gap A, which is a measure of crys- 
tal field splitting, and the Stoner factor /, a measure of 
the Hund coupling. As long as A > / and x\ ^ 1, a 
large region of half-metallicity exists. Both conditions 
are related primarily to the atomic characteristic of con- 
stituents and gross features of the crystal structure, and 
are not sensitive to details. 

The same competition between the band energy and 
the Stoner energy leads to very different magnetic prop- 
erties in case of stoichiometric C0S2, and of C0S2 doped 
with Ni [x > 1 region in Fig.l). The experimental mo- 
ment in the stoichiometric compound is 0.85-0.9 /is, in 
perfect agreement with full potential LSDA calculations, 
and in reasonable agreement with LMTO results as well. 
External pressure, or substituting S by Se, rapidly re- 
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duces magnetic moment ||. As shown in Ref. |l(|, Se 
doping increases the width of the conductivity band (be- 
cause the Se2 dimer has smaller ppa splitting than S2, 
and therefore the Co e q states are better aligned with the 
chalcogen ppa* states |13|). Thus Se doping has the same 
effect as applying pressure. Doping with Ni has similar 
effect and Ni concentrations of the order of 10% make 
the system non-magnetic. The large magnetic moment 
of C0S2 suggests that the ferromagnetism in this system 
is robust and the fact that it disappears so rapidly with 
pressure and doping seems surprising. 




FIG. 2. Density of states for paramagnetic C0S2, calcu- 
lated by LAPW method. E — Ef = 0.35 eV corresponds 
to Coo.9Nio.1S2 in the rigid band approximation. 

Since all these effects are reproduced very well by reg- 
ular LSDA calculations (Fig.l), ESM calculations again 
provide valuable physical insight. Unlike the parabolic 
band case considered above, now one has to take into ac- 
count specific structure of DOS of C0S2 (Fig.|). The loss 
of the one-electron energy can be expressed [ 171 in terms 
of the average DOS, N(M), defined as M/(fj, B H xc ), 
where {^bH xc ) is the (rigid) exchange splitting produc- 
ing magnetic moment M; AE = f Q M MdM/2N(M). 
The best visualization of the ESM is via plotting N(M) 
as a function of M. Wherever this curve crosses the line 
1/7, one has an extremum of the total energy. If the 
slope of the curve is positive at the intersection point, 
the extremum is a maximum, and the state is unstable, 
otherwise it is a minimum and indicates a (meta)stable 
magnetic state. The ESM plot for C0S2, based on the 
LAPW DOS (FigJ|) is shown on Fig.|. At normal pres- 
sure, there are potentially two metastable states: a low 
spin, of the order of 0.3 fiB, and a high spin, of the order 
of 0.85 [Ib- For Stoner parameter I < 0.63 eV neither 
state is (meta)stable, for 0.63 < I < 0.71 eV only the 
low-spin state is, for 0.70 < I < 0.73 eV there are two 
metastable states, and for yet larger I only the high spin 
state may be realized. For I > 0.77 eV ESM produces a 
HM. One can estimate I for a compound using the proce- 
dure described in Ref. pM , or deduce it by comparing the 



ESM to fixed moment calculations. The latter method 
leads to I = 0.68 eV for LMTO and I = 0.72 eV for 
full potential calculations. The complicated structure of 
N(M) can be traced down to the structure of DOS near 
the Fermi level (Fig.|J): after the initial drop of N(M), it 
starts to increase again when the initially fully occupied 
peak at E — Ep w —0.1 eV becomes magnetically polar- 
ized. The high-spin solution corresponds to the situation 
where this peak is fully polarized. 
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FIG. 3. Extended Stoner plot for C0S2, showing the ef- 
fective (averaged) density of states N(Ef) as a function of 
magnetic moment M. The horizontal line corresponds to 
a Stoner factor / = 0.72 eV. The two crossing points near 
M ~ 0.3 fiB, and M « 0.8 fiB, show two (meta)stable 
states in the rigid band approximation. 

Let me now consider the effect of pressure. Roughly 
speaking, applying pressure amounts to rescaling the 
band structure proportionally to squared inverse lattice 
parameter. Correspondingly, the density of states is 
rescaled proportional to (Vq/V) 2 ^ 3 . The total energy in 
the ESM will include that as 




E(V, M) « £ + 
IM 2 



V 
Vo 
B(V 



2/3 



M 



MdM 



2N(M,V ) 
Vo) 2 



2Vr. 



(i 



(1) 



where Vo, Eq, and B are the equilibrium volume, energy, 
and the bulk modulus, respectively, of the paramagnetic 
phase. The magnetic energy shifts the energy minimum 
towards larger volumes, a standard magnetostriction ef- 
fect. What is unusual about the equation of states (|l|) 
with N(M, Vo) from Fig.|| is that for some range of val- 
ues of I there are two local minima, a high spin state 
with a larger volume, and a low spin state with a smaller 
volume. The actual I seems to fall into this range. This 
suggests a first order phase transition with pressure, and 
there are indications p9| that it has been observed in 
the experiment S, although the low spin state, which 
has M w 0.2 ^lb in ESM and « 0.1 /is in fixed mo- 
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ment calculations, was reported in Ref. j|] to have no or 
very small magnetic moment. Another consequence of 
the physical picture outlined here is metamagnetism: In 
the low spin state close to the critical pressure the sys- 
tem can be switched over to the high-spin state by an 
external magnetic field defined by the energy and mag- 
netic moment difference between the two states. Again, 
metamagnetic behavior has been observed in Co(S,Se)2, 
in Coi_ a; Ni :E S2, and in compressed C0S2. Finally, the the- 
ory predicts rapid increase of the equilibrium magnetiza- 
tion with negative pressure: in the ESM an expansion of 
3-4% in volume already increases M to nearly 1 hb- This 
suggests that the spin fluctuations at high temperature 
(above Tc) may have larger amplitude than the ordered 
moments at zero temperature. This kind of behavior has 
also been observed ^0|, and discussed in the literature 

ID- 

The family of pyrite materials formed by 3d transition 
metals and chalcogens is incredibly rich. It shows various 
kinds of magnetism and metamagnetism, metal-insulator 
transitions of different types, superconductivity, and half 
metallicity. Except for the vicinity of a Mott-Hubbard 
transition, that is, close to N1S2, the physics of these ma- 
terials can be rather well understood within the local spin 
density functional theory. In particular, the extended 
Stoner formalism provides considerable insight into the 
magnetic behavior of this system. The most important 
conclusions from the calculations are: 

1. The Fei-zCoa^ alloy is predicted to be a half 
metal in a large range of concentration. Unlike most 
other known half metals, and very importantly for the 
applications, the half-metallicity should be robust with 
respect to defects and crystallographic disorder. 

2. C0S2 is an itinerant fcrromagnet which has, due to 
a complicated structure of its density of states, two mag- 
netic states: a high-spin one, with the moment M ~ 0.8 
jj,B, and a low-spin one, with M < 0.1 [Ib- The first order 
transition to the low spin state can be induced by exter- 
nal pressure, doping with Se, or with Ni. In the low spin 
state the material exhibits metamagnetic properties. The 
structure of the density of states also manifests itself via 
an unusual temperature dependence of the magnitude of 
magnetic fluctuations. 

3. The closeness of C0S2 to a magnetic phase transition 
also leads to strong coupling between electronic, lattice, 
and magnetic degrees of freedom, making it a relative of 
such "bad metals" as SrRu03 and magnetoresistive man- 
ganites. Therefore interesting transport properties are 
to be expected in the system at and near stoichiomet- 
ric composition, including large magnetoresistance and 
violation of the Yoffe-Regel limit at high temperatures. 

For such an interesting, the 3d pyrites seem to be un- 
usually little studied either experimentally or theoreti- 
cally. Hopefully this paper will encourage further inves- 
tigations. 
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